Design and Stability of Discrete-Time Quantum 
Filters with Measurement Imperfections 

Abhinav Somaraju Igor Dotsenko Clement Sayrin Pierre Rouchon 



Abstract — This work considers the theory underlying a 
discrete-time quantum filter recently used in a quantum feed- 
back experiment. It proves that this After taking into account 
decoherence and measurement errors is optimaf and stable. We 
present the generaf framework underlying this After and show 
that it corresponds to a recursive expression of the feast-square 
optimal estimation of the density operator in the presence of 
measurement imperfections. By measurement imperfections, we 
mean in a very general sense unread measurement performed 
by the environment (decoherence) and active measurement per- 
formed by non-ideaf detectors. However, we assume to know 
preciseiy all the Kraus operators and also the detection error 
rates. Such recursive expressions combine well known methods 
from quantum Altering theory and classical probability theory 
(Bayes' law). We then demonstrate that such a recursive Alter is 
always stable with respect to its initial condition: the Adelity 
between the optimal Alter state (when the initial Alter state 
coincides with the real quantum state) and the Alter state (when 
the initial Alter state is arbitrary) is a sub-martingale. 

I. Introduction 

The theory of filtering considers the estimation of the system 
state from noisy and/or partial observations (see, e.g., [5]). For 
quantum systems, filtering theory was initiated in the 1980s by 
Belavkin in a series of papers [1], [2], [3], [4] (also see the tu- 
torial papers [7], [6] for a more recent introduction). Belavkin 
makes use of the operational formalisms of Davies [9], which 
is a precursor to the theory of quantum filtering. He has also 
realized that due to the unavoidable back-action of quantum 
measurements, the theory of filtering plays a fundamental role 
in quantum feedback control (see e.g. [4], [2]). The theory of 
quantum filtering was independently developed in the physics 
community, particularly in the context of quantum optics, 
under the name of Quantum Measurement Theory [8], [12], 
[11], [16]. 

Most of this theory has been developed for continuous-time 
systems and little emphasis has been given to measurement 
imprecisions and their explicit impact on the filter design 
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and time-recursive equations. To our knowledge, the problem 
of designing a quantum filter in the presence of classical 
measurement imperfections has not been examined in the 
discrete time setting. In this paper, we focus on this issue and 
propose a systematic method to derive quantum filters taking 
into account several detection error rates. 

In [11, Sec. 2.2.2], the authors discuss how the state of 
a quantum system evolves after a single imprecise measure- 
ment. In [10], a recursive quantum state estimation with 
measurement imperfections has been considered. In [14], such 
a quantum state estimate has been used in a quantum feedback 
experiment that stabilizes photon-number states of a quantized 
field mode, trapped in a super-conducting cavity. We prove 
here that such estimates are in fact optimal since they coincide 
with the conditional expectation of the quantum state (density 
matrix) knowing the past detections, the error rates and the 
initial quantum state. 

Section II describes the structure of a genuine quantum 
measurement model including detection error rates which 
is a straightforward generalization of the models considered 
in [11], [10], [14]. This model may be used in situations 
with partial knowledge of all the quantum jumps and also 
measurement errors of the jumps that are detected. However, 
they assume to know precisely all the Kraus operators and also 
the detection error rates. 

Section III is devoted to the first result in this paper 
summarized in Theorem 111. 1 : the conditional expectation of 
the quantum state knowing the past detections and the initial 
state obeys a recursive equation in each discrete time-step. 
This recursive equation is given in (3) and depends explicitly 
on the error rates. The proof of Theorem III. 1 shows that such 
recursive equation may be derived by a simple application of 
Bayes' law. 

In section IV, we prove that the quantum filter defined in 
Theorem III. 1 is stable versus its initial conditions: the fidelity 
between the optimal estimate conditioned on the initial state 
of the system being known and a second estimate in which the 
initial state is unknown is a sub-martingale. This stability result 
combines Theorem III. 1 and [13]. Note that stability does not 
imply convergence, in general. For convergence results in the 
continuous-time case see, e.g., [15] and the references therein. 

In section V, we describe in detail the Kraus operators 
and error rates modeling the discrete-time quantum system 
considered in the quantum feedback experiment [14]: the 
quantum filter used in the feedback loop corresponds precisely 
to the recursive equation (3) given by Theorem III. 1; according 
to Theorem IV. 1, this filter tends to forget its initial condition. 
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II. Measurement Model 

In this section we discuss the model describing repeated and 
imperfect measurements on a quantum system. Such modeling 
including decoherence-induced quantum jumps and measure- 
ment errors is a direct generalization of the one proposed 
in [10] and used in real-time for the quantum feedback exper- 
iments reported in [14] (also see [11]). We initially consider 
the case of a single ideal measurement and then develop the 
model to consider imperfect and repeated measurements. The 
final model is described in Subsection II-C. 

A. Ideal Case 

Let H be the system's Hilbert space with p\ a density matrix 
denoting the initial state of the system at step k = 1. We 
consider the evolution p\ i-> p2 of such a quantum system 
under discrete-time quantum jumps (see e.g. [12, Ch. 4] or [11, 
Ch. 2]). 

Consider a set of Kraus operators M q : H — > H, q £ 

M^M q = I. Here we assume 
£ N possible quantum jumps and I is 
the identity operator on H. The superscript id stands for an 
abbreviation of ideal. Consider an ideal world with full access 
to all quantum jumps via a complete and ideal set of jump 
detectors. When the quantum jump indexed by q is detected, 
the state of the system changes to 
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is the probability to detect jump q, knowing the state p\. We 
now consider the case of realistic experiments with possible 
measurement errors. 

B. Realistic Case (with imprecise measurements) 

We consider that the ideal detection of the jump q corre- 
sponds to an ideal measure outcome p ld = q. We denote 



by random variable p ld £ {1,2, 



i } this outcome 



provided by ideal sensors. We assume that realistic sensors 
provide an outcome /i rl that is a random variable in the set 

{1,2,..., 777^}. 

We assume that, with a known probability, an ideal mea- 
surement outcome /i ld occurs effectively whereas the realistic 
sensors detect an outcome pf 1 . The correlations between the 
events p ld = q and p rl = p are modeled by classical 
probabilities through a stochastic matrix 77 £ IR m xm : 

Vp , q =V[p* 1 =p\p ld = q}. 

It gives the probability that the real sensors detect p zl = p 
given the ideal sensors would detect p ld = q, for p £ 
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<7> E p =i Vp,q = 1> tne matrix 77 = (r) Pt q) is a left stochastic 
matrix. 

C. Realistic Experiment with Repeated Measurements 

Consider the case of a sequence of discrete-time mea- 
surements. We denote by p k the state of the system at 
discrete time-step k. Also suppose M q - k is the Kraus op- 
erator corresponding to the k th ideal measurement for q £ 
{1, 2, . . . , m ld }. Note that we allow for a different set of 
Kraus operators Al.-k for different time-steps k. One can also 
consider m ld and m rl be dependent on k. 

Similar to the previous subsection, we denote by pSj} £ 
{1, . . . , m rl } and pif? £ {1, . . . , m ld }, the random variables 
corresponding to the k th realistic and ideal outcomes, respec- 
tively. Therefore, 
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andP[/4 d = <7 I Pk] =Tr- 

Also, we assume rf £ IR™***™"*, the stochastic matrix de- 
termining the probability of error, can depend on the discrete- 
time step k. In particular, we have 

n^t=pH i = q] = Vp, q , 

for p £ {1, . . . , m rl }, q £ {1, . . . , m ld }. 

III. Recursive Equation for the optimal filter 

We wish to obtain a recursive equation for the optimal 
estimate f>k+i of the state p k +i knowing initial value p\ 
and real measurement outcomes pT^, 
estimate pk is defined as 



/n^ 1 . This optimal 



pk =E\p k \p 1 ,pf-,.. 
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The following theorem says that we can ignore the original 
state pk and only consider p k that is shown to be the state of 
a Markov process. 

Theorem III.l. The optimal estimate pk satisfies the following 
recursive equation 



Pk+i = 



E7^vt q M q;kPk Mi k 



ifl€ 



Pk- Moreover, we have 

Pl,---,p T k -i =Pk-l 



(3) 
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Remark III.l. The division in (3) by the R.H.S of (4) 

could appear problematic when this denominator vanishes. 



Nevertheless, if we assume that the real measurements are 
Pi, 



pf 



i pf = Vk, then 
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(otherwise such measurement outcomes are not possible). 
Consequently, 
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cannot vanish. Thus recurrence (3) is always well defined 
because we have (4). 

Remark III.2 (Markov property of the filter). Equations (3) 
and (4) tell us that the joint-process (p k , Pk ) is a Markov 
process and therefore the statistics of the measurement process 
pf may be determined using p k . This in particular implies 
that we may use Monte Carlo methods to simulate the obser- 
vation process pf only using p k independent of the actual 
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state pk and measurement history pf , . . 



Proof: In this proof we use the following notation for 
ease of presentation: we use pf = q L and pf = p L to 
denote the set of events {p\ d = qi,pf = qi, . . . , pf = q k } 
and {pf = pi,pf =P2,---, pf — Pk}, respectively. For in- 
stance, using this notation, we have 
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Assume that the values measured by the real detector are 
Pi = pf, pf. = pf. Then we have the optimal estimate 



Pk + 1 
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Consequently, we have 

p [pf = 1l\pi, pf = Pl] 
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Injecting the above relation into (5) and using (6) yields 

Pfc+i = 

ggi qk <^i ■ ■ ■ ri k . qk M *k* ■ ■ ■ M, i; ipiA/t i:1 ■ ■ ■ Ml h . k 

(9) 

It is then clear that p k +\ can be calculated from p\ in a 
recursive manner according to 

T. q V k Pk , q M q , kPk Ml. k 

pk+1 H^k-^'-^KA' 



where the intermediate states correspond to optimal estimates 
from step 2 to b. p 3 = E[ P j\ P i, pf = pi, . . ■ , pf-i = Pj-i], 
j = 2, . . . , k. The recursive relation (3) is thus proved. 

We now prove (4). In the following, we set ordered product 

M, 4 M qj;j ■ ■ ■ M q2 . 2 ■ M qul 



for j = l,...,fc. 
We have 
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Pi,pf = Pi, • ■ • , pf-i =Pk-i 
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This fraction can be computed using (8): 
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According to (9) with k — 1 instead of k, we have 
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we get finally (4). 



IV. Stability with respect to initial conditions 

Assume that we do not have access to the real initial state 
p\. We cannot compute the optimal estimate pk- We can still 
use the recurrence formula (3) based on the real measurement 
outcomes (/ij 1 = Pj)j=i,..., k -i to propose an estimation p k 
of pk- We will prove below that this estimation procedure is 
stable in the sense that the fidelity between p k and p k is non- 
decreasing in average whatever the initial condition p\ is. 

For ease of notation we set 



M 



p,q;k 



M p . k (p) 



q;k, 
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for any k > 1, p £ {1, . . . , m rl } and q E {1, . 
sets 

Sp-k — {MpS;k, ■ ■ ■ , M p m ±d. k } 

for p = 1 , . . . , m rl form a partition of 

Sk = {Mp^. k | p = 1, . . . , m rl , ? = 1, . 
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i }. The 
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Using this notation, the recursive Equation (3) defines a 
coarse-grained Markov chain in the sense of [13]. 

If p r k = p, we define p\ recursively for k > 1 as follows 
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a recursive formula is valid as soon as 
Eg'Tr { M p,q;kPlMp,g;k} > °> wnicn is automatically 
satisfied when p| is of full rank. p% +1 is indeterminate when 

Ttql* { M P,<r,kP% M l,q;k} = °- But usin g the continuity 
arguments developed at the end of the appendix we can give 
a value for p% +1 in the following way: for each p k , consider 
the set of density operators 

KtK £ )) H ,..., mI11 

where e > and ^ e = + eI)/Tr {p^ + el}. Since 
p% e is positive definite, ( M. v -u (p% e ) ) are well 

\ ' / p=l,...,m rl 

defined and admit a limit point when e tends to + (W is of 
finite dimension here). Take for each p k such a limit point 

(Pk+l, p )p=h-,m*- Set Pfc+1 = Pt+hP When Mfc = P- If 



J2 q Tv { M P,q;kPk M lq-k\ > then Pfc+i.p coincides neces- 
sarily with -Mp-k(pl.) and we recover (11). During the proof 
of Theorem IV. 1, we will use only recurrence (11) having in 
mind that, when E g r ^ r \j^p,i;kPk^p,q;k\ = we have to 

use Pt+i = Pt+i, P - 

If the initial state of the system p\ coincides with pi, 
then pf. coincides with the optimal estimate p k of the state 
p k from Theorem III. 1 . In fact, once the initial states p\ 
and p\ are given, the process (p k ) and (pi) are driven 
by the same stochastic process (p^-^), itself driven by the 
combination of the original quantum process of state (pk) with 
the classical process associated to the left stochastic matrices 
(rj k ) governing detection errors. 

The following theorem shows that the fidelity between p k 
and p\ is non-decreasing in average. 

Theorem IV. 1. Suppose p e k satisfies the recursive relation (11) 
with an arbitrary initial density operator p\. Then the fidelity 
between p k and pt defined by 



'PkPkVPk 



is a submartingale in the following sense: 
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Proof: Denote Fk = F(pk,P k )- We have 
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The conditional probabilities appearing in this sum are given 
by (4) and the conditional expectations read 



F k 
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since, once p\ and p\ and pJ^ = p\,..., p r k-i = Pk-i and 
pi 1 = p are given, p k+1 = M p -k(pk) and P k +i = ^pApD- 
Thus we have 
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direct consequence of equation (12) given in appendix with 



some operators that are used to describe these quantum jumps: 

D a = e aaf - a * a , L no = V-Pa(O) I, L g = VPa(l)cos0 N , 



id , index j corresponding to p, index i Lp = ^pj^j sin ^ Lgg = ^pj^ cog 2 ^ 



to (p,q), operators Li to M p ^ q: k, density operators p and a to 
p/c and /o|, respectively. 



V. Example: Quantum Filter for the Photon-Box 



This section considers, as a key illustration, the quantum 
filter design in [14] to estimate in real-time the state p of 
a quantized electro-magnetic field. Since this filter admits 
exactly the recursive form of Theorem III. 1 , Theorem IV. 1 
applies and thus, this filter is proved here to be stable versus 
its initial condition. 

The actual experiment under consideration uses quantum 
non-demolition measurements [12] to estimate the state of the 
quantized field trapped in a superconducting microwave cavity. 
Circular Rydberg atoms are sent at discrete time intervals to 
perform partial measurements of the photon number. Atoms 
are subsequently detected either in their excited (e) or ground 
(g) state. The outcomes of these measurements are then used 
to estimate the state of the cavity field, thanks to the quantum 
filter described below. This estimation is eventually used to 
calculate the amplitude of classical fields injected in the cavity 
in order to bring the field closer to a predefined target state. 
The interested reader is directed to [10] and [14] for further 
details of the experimental setup and results obtained. 

The Hilbert space T~L of the cavity is, up to some finite 
photon number truncation (around 10), the Fock space with 
basis {|n)} n >o, each \n) being the Fock state associated to 
exactly n photons (photon-number state). The annihilation 
operator a : % n- % is defined by a |n) = \fn\n — 1) for 
n > 1 and a |0) = 0. Its Hermitian conjugate a^ is the creation 
operator satisfying a |n) = \/n + 1 \n + 1), for all n > 0. The 
photon-number operator (energy operator) is N = a^a which 
satisfies N |n) = n \n), for all n > 0. Recall the commutation 
[a,at]=I. 

In [14] the following imperfections have been considered: 

• atomic preparation efficiency characterized by P a (n a ) > 
0, the probability to have n a € {0, 1, 2} atom(s) interact 
with the cavity: P a (0) + P a (l) + P a (2) = 1. 

• detection efficiency characterized by € [0,1], the 
probability that the detector detects an atom when it is 
present. 

• state detection error rate r\ g € [0, 1] (resp. ?y e ) probability 
of erroneous state assignation to e (resp. g) when the 
atom collapses in g (resp. e). 

The original state p is subject to m ld = 3x7 possible 
quantum jumps and the available sensors (atomic detector) 
admits only m rl = 6 possibilities. We begin by introducing 



L ge = L eg = \/Pj2) cos<?!) N sin0 N , L ee = V /P a (2)sin 2 </> N , 



L n = l- 



t(l+2n th ) 



N 



where 



N = 0o(N+1 2 /2)+0fl and < e,n th « 1, o ,0 fl are 
real experimental parameters. The unitary displacement oper- 
ator D a corresponds to the control input a E C, depending on 
the sampling step k. The operators L a , L + and L_ correspond 
to the interaction of the cavity-field with its environment 
(decoherence due to mirrors and thermal photons): 

1) L corresponds to no photon jump; 

2) L + corresponds to the capture of one thermal photon by 
the cavity-field; 

3) L_ corresponds to one photon lost from the cavity-field. 
Since L\L + L\_L+ + L*_L- = I + 0(e 2 ) and e is small, we 
consider in the sequel that (L Q , L + , L_) are associated to an 
effective Kraus map L ph\ + L + pL* + + L-pl}_. 

The operators L no , L g , L e , L gg , L ge , L eg and L ee corre- 
spond to the jump induced by the collapse of possible crossing 
atom(s) having interacted with the cavity-field: 

1) L no - no atom in the atomic sample; 

2) L g - one atom having interacted with the cavity-field and 
collapsed to the atomic ground state during the detection 
process; 

3) L e - one atom having interacted with the cavity-field and 
collapsed to the atomic excited state during the detection 
process; 

4) L gg - two atoms having interacted with the cavity-field, 
both having collapsed to g; 

5) L ge - two atoms having interacted with the cavity-field, 
the first one having collapsed to g and the second to e; 

6) L eg - two atoms having interacted with the cavity-field, 
the first one having collapsed to e and the second to g; 

7) L ee - two atoms having interacted with the cavity-field, 
both having collapsed to e. 

For each control input a, we have a total of ?n ld = 3x7 
Kraus operators. The jumps are labeled by q = {q a ,q c ) with 
q a G {no,g,e,gg,ge,eg,ee} labeling atom related jumps 
and q c G {o, +, — } cavity decoherence jumps. The Kraus 
operators associated to such q are M q = L q cD a L q a. So, for 
instance, with the control input au at step k: 

• the Kraus operator corresponding to one atom collapsing 
in ground state, q a = g, and one photon destroyed by 
mirrors, q c = — , reads M^ g _■. = L-D ak L g . 

• the Kraus operator corresponding to two atoms, the first 
one collapsing to g, the second one to e, q a = ge, and one 
thermal photon being caught between the two mirrors, 
q c = +, reads M^ ge +) = L + D ak L ge . 

One can check that, for any value of a, these 21 operators 
define a Kraus map (using the assumption that L\L +L\L + + 
LlL- w I). 
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TABLE I 
Stochastic matrix 77 Pj9 . 



We have only m 11 = 6 real detection possibilities p <G 
{no,g,e,gg,ge,ee} corresponding respectively to no detec- 
tion, a single detection in g, a single detection in e, a double 
detection both in g, a double detection one in g and the other 
in e, and a double detection both in e. A double detection does 
not distinguish two atoms. The entries r) Ptq of the stochastic 
matrix describing the imperfect detections corrupted by errors 
are independent of a and given by Table I. It relies on error 
rate parameters r\ g , rj e and ed in [0, 1] and assume no error 
correlation between atoms. So for instance the probability that 
there is a single atom in the sample, which collapses to the 
ground state and is in fact detected by the experimental sensors 
to be in the ground state is given by r)g,( g ,qc) = £d(l — 
Note that each column in the table sums up to one, since rj P} q 
is a left stochastic matrix. 

VI. Conclusion 

In this paper, we derive a recursive expression for the 
optimal estimate of a quantum system's state from imperfect, 
discrete-time measurements. The optimality of this recursive 
expression is proven by a simple application of Bayes' lemma 
and quantum measurement postulates. Such a filter is shown 
to satisfy a Markov property and thus can be used for quantum 
control as shown in [14] or to run Monte Carlo simulations 
of the measurement trajectories. We also demonstrate that this 
filter is stable with respect to initial conditions in the sense 
of Theorem IV. 1 . In the future, the continuous-time version of 
these results will be investigated. 
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Appendix 

Consider a set of s operators (ij)i=i,..., s on the finite 
dimensional Hilbert space H, such that Li ^ for all i and 
Y^i—i L[Li = I. Take a partition of {1, ...,s} into r > 1 
non-empty sub-sets {Vj)^_ x r - Then, for any semi-definite 
positive operators p and a on H with unit traces, the following 
inequality proved in [13] holds true: 

r 

F(p,a) < ^Tr{E ie7=j .L lP L]}x 

< { : : ■ ; " ; . ' " d2) 

^{Z^.LipLl} Tr{-Z ieVj L it ,L\} J' 

where F(p, a) = F(cr, p) is the fidelity between a and p 
defined as 

F(p,cr) = (TrjyV^V^}) • 

When, in the above sum, a denominator Tr jX^ep 
depending on p vanishes, the sum remains still valid since the 
corresponding value of F is multiplied by the same vanishing 



factor and F is bounded since between and 1. This is no 
more true when a denominator Tr ^i "^!} depending 

on a vanishes. In this case the above inequality should be 
interpreted in the following way relying on a continuity 
argument. For each a, define 



Z a = 




Li<TL\\ =0) C {!,... ,r}. 



For almost all density operators a, Z a — 0. Take e > and 
consider the positive definite density operator er e = ^r^i^n ■ 

For j G {1, • ■ • ,r}, Tr {E;e^ L i^ L l} > since each L< ^ 
0. In particular, the set of density operators 



admits at least an accumulation point when e tends to + (H 
of finite dimension implies that the set of density operators is 
compact). Denote by {o't)j£Z <T , such an accumulation point 
where each aj is a density operator. Since for any e > 0, 
inequality (12) holds true for p and a e , F is continuous, and 

* ■ A -3 L t a'L\ £ UuL\ 

tor any i e Z„, — ; i r-r tends to — , - tt, 

we have by continuity 

Such continuity argument show that we can extend inequal- 
ity (12) via the accumulation value(s) er+ when Z a is not 
empty. 



